Effect of Positional Disorder on Transmissional Properties of a Two-Dimensional 

Photonic Crystal: An FDTD Study 



A. R. Hashemi, M. Hosseini Farzad, and Afshin Montakhab* 
Department of Physics, College of Sciences, Shiraz University, Shiraz 71454, Iran 

(Dated: September 7, 2009) 

In this paper we study the effect of positional randomness on transmissional properties of a 
two dimensional photonic crystal as a function of a randomness parameter a {a — completely 
ordered, a = 1 completely disordered). We use finite-difference time-domain (FDTD) method to 
solve the Maxwell's equations in such a medium numerically. We consider two situations: first a 
90° bent photonic crystal wave-guide and second a centrally pulsed photonic crystal micro-cavity. 
We plot various figures for each case which characterize the effect of randomness quantitatively. 
More specifically, in the wave-guide situation, we show that the general shape of the normalized 
total output energy is a Gaussian function of randomness with wavelength-dependent width. For 
centrally pulsed photonic crystal, the output energy curves display extremum behavior both as 
a function of time as well as randomness. We explain these effects in terms of two distinct but 
simultaneous effects which emerge with increasing randomness, namely the creation of semi-localized 
modes and the shrinking (and eventual destruction) of the photonic band-gaps. Finally, we attribute 
these effects to the overlap probability of random dielectric cylinders within the photonic crystal as 
randomness is increased. 

PACS numbers: 42.25.Dd, 42.70.Qs, 78.40.Pg, 78.20.Bh 



I. INTRODUCTION 

Photonic crystals (PC), structures with a periodic re- 
fractive index distribution, are a subject of intense world- 
wide research In contrast with the traditional wave- 
guides in which guiding principle is based on total in- 
ternal reflection, PC wave-guides have a fundamentally 
different mechanism''. These wave-guides can propagate 
the light in a path on the scale of wavelength without 
diffraction losses even at corners, because light leakage 
is forbidden for frequencies in the photonic band-gap^. 
Among various PC structures, two dimensional (2D) PCs 
based on a finite cluster of cylinders of infinite length 
have attracted most attention-''^, since these structures 
are relatively easy to fabricate or simulate. 

The effect of disorder in 2D PCs with 90° bent wave- 
guide was investigated by Langtry et al. ^ and also by 
Kwan et al. ^. In the former, the sensitivity of wave- 
guiding to different degree of disorder mostly in radius of 
cylinders, refractive index, and filling factor have been 
studied by highly accurate multipole method. In the 
latter, by using multiple-scattering method, it has been 
shown that the transmission spectrum of the wave-guide 
is more sensitive to the uniformity in position of first 
layer surrounding the wave-guide. 

In the first part of this article, we study, using finite- 
difference time-domain (FDTD) method*, normalized 
output energy of a 90° bent wave-guide in 2D PC consist- 
ing of a square array of long dielectric cylinders in the air 
background. By introducing a wide range of randomness 
in the position of all cylinders, we show that the output 
energy normalized to the ordered case have a Gaussian 
shape, when plotted with respect to randomness factor. 
This behavior is shown for any wavelength in the original 
PC band-gap. 



In the second part of this article, we study the effect of 
positional randomness on the light energy output from a 
PC micro-cavity, where four cylinders are removed from 
the center of square array and a short Ricker wavelet 
pulse® is initiated in the center of the micro-cavity. We 
use standard FDTD method to study the time evolu- 
tion of output energy for TM electromagnetic field in 
the above system. The normalized output energy curve 
as a function of time for different randomness factors 
starts out with a sharp minimum followed by a broad 
maximum. Previous investigations are mainly concerned 
with cavity mode frequencies and quality factor versus 
strength of randomness^^*^ and on the variation of density- 
of-states with respect to reduced frequency^-'. Also, some 
researchers have undertaken theoretical studies of po- 
sitional disorder effect on resonant mode distributions 
and resonance frequencies in 2D photonic crystal micro- 
cavity'^''^, and experimental activity accompanied by 
numerical simulation of spatial intensity distribution of 
localized modes in 2D open microwave cavity'''. In the 
latter, for detection and simulation of localized modes, 
the random medium is excited by more than one thin 
antennas around the center of randomized PC. However, 
in the present work, we excite the medium by only one 
central Ricker wavelet pulse and then the system is left 
to its own in order to study how the energy of the pulse 
exit from boundaries of PC when it is subjected to the 
different strengths of randomness. 

Using transfer matrix method, the first study of the ef- 
fect of induced randomness in radius or refractive index of 
the cylinders have shown the appearance of propagation 
states in the original band-gap'^. Band-gaps, although 
with some narrowing, were found to persist even for large 
amount of disorder. This phenomenon has been observed 
both theoretically'^"'® as well as experimentally^°'^', es- 
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pecially in the case of positional disorder in cylinders^^'^^. 
These investigations mainly concerned analysis of trans- 
mission through PC slabs. However, in the present work 
the source is placed in the center of our system and band- 
gap resistance to disorder is observed in the transmission 
spectrum of the output light. We show that the forma- 
tion of non-central semi-localized modes, besides the light 
localization in the central micro-cavity, can explain dy- 
namical behavior of normalized output energy for various 
randomness factors. 

This paper is arranged as follows: in Section II wc 
briefly explain the set-up for our numerical simulations. 
Section III entails the main results of our study which 
includes numerical results as well as discussion of the 
relevant physical effects. We summarize our main results 
in Section IV. 



domness. For first scenario, i.e. propagation of light by 
a wave-guide in PC [Fig. 1(c)], we enter a plane- wave to 
the wave-guide and calculate the output energy at the 
other end; while for the second scenario, i.e. propaga- 
tion of a pulse from the center of PC, we use a Ricker 
wavelet pulse with central wavelength, A — 550 nm, as 
initial pulse and calculate total output energy on all sides 
of crystal. In the second scenario, we remove four cen- 
tral cylinders to prevent the source from overlapping with 
randomly arranged neighboring cylinders and we there- 
fore have a cavity in the center of the PC [Fig. 1(b)]. 

III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. wave-guide PC 



II. NUMERICAL SIMULATION 

In this paper we concentrate on a PC model which 
consist of a square lattice of cylinders with refractive in- 
dex 3 with radius 0.3d, where d is the lattice constant 
of the crystal [Fig. 1(a)]. The surrounding of cylinders 
is air with refractive index 1. This kind of crystal has a 
band-gap in wavelength region 3.0(i < A < 3.7(i^^. We 
consider the cylinders axes aligned in z direction, so in 
order to introduce disorder in the position of cylinders, 
we impose displacements bdx and 8dy for each cylinder in 
X and y directions respectively, with the random magni- 
tudes between zero and a percentage of lattice constant 
d. We call this percentage randomness factor a. Thus in 
the disordered case, like Fig. 1(b), each cylinder displaces 
from its ordered position, by Sdx = Ux x a x d in x direc- 
tion and by Sdy = Uy x a x d in y direction; where Ux and 
Uy are random numbers with flat distribution between -1 
and +1. We perform our simulation with d = 180 nm 
and so, we choose visible light, A ^ 550 nm, in the edge 
of band-gap, to study the effect of randomness better. 
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FIG. 1: (a) Perfect photonic crystal, (b) disordered PC {a = 
0.5) with cavity, (c) disordered PC (a — 0.5) with wave-guide. 



Figure 2 shows normalized output of guided energy in a 
90° bent wave-guide PC. For each randomness factor, the 
output energy is averaged over 200 different random sam- 
ples and normalized to the output energy of the ordered 
case (a — 0). As we expect, output energy decreases 
with increasing disorder, however, we note the following 
observations: 

i) The general shape of the output energy curve 
seems to fit well with a Gaussian function, i?oMt = 

2/2 

EorderG"" ■ Increasing number of samples used 
in averaging makes an increasingly better fit with 
such function. 

ii) We observe that the general Gaussian shape is in- 
dependent of the wavelength of light we used [Fig. 2 
insets]. However, the width of the Gaussian (cr) de- 
creases with increasing A, whether A lies in the edge 
of band-gap or not. 

in) The relatively large value of a indicates certain ro- 
bustness to disorder in this wave-guiding effect. For 
example, in Fig. 2, at 15% randomness, output en- 
ergy is still about 90% of the ordered amount. 

It is obvious that the above observations is related to 
the relative robustness of the photonic band-gap struc- 
ture for small and intermediate values of a and its even- 
tual destruction for large a [see Fig. 7]. However, the 
Gaussian shape of the figure is a curious effect, a point 
which we will come to later on this paper. 



To simulate propagation of light, we use FDTD 
method, and for truncating computational space, we use 
CPML absorbing boundary condition^^. 

Our model is symmetric along z axis, so choosing TMz 
mode, we carry out the simulation in two dimensions. For 
each randomness factor a, we repeat simulation with sev- 
eral different random arrangements (ensemble) of cylin- 
ders, and average over the various realizations of ran- 



B. centrally pulsed PC 

We now turn our attention to the second scenario 
where a Ricker wavelet source is placed in the middle 
of the PC and the total output is monitored along all 
edges. To study the changes in localized modes in the 
crystal, we plot output energy of pulse from the crystal 
as a function of time for different randomness factors in 
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FIG. 2: Output energy of guided light by a 90° bent wave- 
guide in the photonic crystal (solid line) compared with a 
Gaussian function with a — 0.392 (dashed line). Upper inset 
shows the output for A — 570 nm compared with a Gaussian 
with a = 0.327, and lower inset shows the output for A = 600 
nm compared with a Gaussian with a — 0.262. The axis on 
the inset are the same as the main figure. 

Fig. 3. For each a, we average output energy over 300 
different random samples and subtract output curve of 
ordered crystal (a = 0) from it, then normalize the re- 
sult to the total input energy. Therefore, in Fig. 3, each 
curve represents the (average) effect of a given random- 
ness (a) on the relative output energy as a function of 
time. An important characteristic stands out: each curve 
starts out with a sharp minimum followed by a broad 
maximum before setting into a constant zero value. 




FIG. 3: Normalized output energy of crystal as a function of 
time for different randomness factors. The output energy of 
the ordered crystal has been subtracted. 

As is shown in Fig. 4, the (absolute) value of the min- 
imum increases with increasing a, while the peak of the 
broad maximum is largely unaffected by the random- 
ness. On the other hand, the minimum occurs at a fixed 
time regardless of a, while the maximum occurs at longer 



times with increasing a (inset Fig. 4). 
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FIG. 4: Energy of maximum and minimum of Fig. 3 ver- 
sus randomness factor. Inset shows time of occurrence of the 
maximum and minimum. Note that the maximum is insensi- 
tive to a while the minimum decreases continuously with a. 
On the other hand, the time of occurrence of the minimum 
is insensitive to a while the peak of the output shifts with 
increasing disorder. 

It is easy to understand the behavior in Fig. 4 by 
considering the effect of randomness on the PC. The 
minimum corresponds to the leaving of the unaffected 
modes from the crystal, i.e. the modes outside the 
band-gap thus leaving the crystal straight away without 
impedance. Therefore the constant time of occurrence 
for this minimum can easily be understood. 

As randomness increases the band-gap shrinks, so we 
expect that a larger portion of the modes leave the crystal 
without being affected; but we see that the minimum 
is lowered by randomness, which means, the unaffected 
modes decreases by increasing randomness. The root of 
this apparent discrepancy is in the creation of the so- 
called semi-localized modes. 

As can be seen from Fig. 5, increasing randomness 
raises the possibility of creation of semi-localized modes. 
Here, in addition to the central cavity localized mode, 
some other modes are confined and remain active inside 
the PC (see the V shape pattern pulsating on the right 
side of the center in Fig. 5). This effect is due to the re- 
peated internal reflections inside the randomly arranged 
cylinders which are somehow coupled or synchronized'^^ . 
Since reflections are not total, these semi-localized modes 
weaken as time elapses. Note that increasing randomness 
allows for more suitable situations to be created for semi- 
localized modes, and increasing the number of these pos- 
sible modes broadens the output curve as seen in Fig. 3 
and also cause the maximum to occur later. Also, it 
is now easy to understand the lowering of the minima 
in Fig. 4: creation of semi-localized modes temporarily 
traps certain modes thus lowering the output energy in 
short times, releasing it later over a broad range of time. 

How is the total energy output affected by increasing 
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FIG. 5: Snapshots of the field in a sample with a = 0.5 at 
different times (blue for positive and red for negative values) . 
Semi-localized mode magnified below. 
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FIG. 6: Total output energy from the crystal, normalized to 
the total output of ordered case, versus randomness factor. 




X/d 



FIG. 7: Transmission of PC for three different randomness 
factors as a function of normalized wavelength (A/d). 



randomness? The answer is shown in Fig. 6 where in- 
creasing randomness increases the total output energy. 
As pointed out above, this is due to the shrinking of the 
photonic band-gap and the resulting effect which allows 
for more modes to escape the PC. The shrinking of the 
photonic band-gap as a result of disorder can be best 

seen in Fig. 7 where the transmission, T — ^ ^°"tp"^* ^ ^ 

has been plotted as a function of wavelength (X/d) for 
three different randomness factors. We have actually 
plotted such figures for various values of a and note that 
the general structure of the primary (wider) band-gap re- 
mains intact up to q; = 0.5 and steadily shrinks for larger 
values of a. In the secondary (narrower) band-gap the 
observed effect of randomness is stronger, however the 
above-mentioned general scheme still holds. 

Creation of semi-localized modes is a direct result of 
coupling of scatterers, so by calculating the probability of 
coupling of cylinders we can better answer the question 



of how randomness affects the propagation of electromag- 
netic waves. In our model, coupling happens when two 
cylinders are placed such that their center to center dis- 
tance is less than or equal to the diameter (2r). We 
therefore calculate the probability of coupling of a cylin- 
der with its eight nearest neighbor in a crystal for a given 
a and average over different realizations of the crystal 
(ensemble averaged). The result is plotted in Fig. 8. 

The shape of this curve is very interesting and we be- 
lieve underlies many of the observed behaviors in our 
study. Two simultaneous but distinct phenomena occur 
as randomness is increased^^: First is the shrinking and 
eventual destruction of the band-gap structure. Second 
is the creation of semi-localized modes. The correlation 
between the general shape of Fig. 8 and many of our 
results (e.g. Figs. 2, 4, 6) is indicative of this general be- 
havior. In Fig. 8, as a is increased from zero no significant 
change occurs, as a sharp (exponential) increase in the 
probability follows in the intermediate a range followed 
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IV. CONCLUSION 
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FIG. 8: Overlap probability of a cylinder with its eight neigh- 
bors versus randomness factor. The solid line represents an 
exponential fit with the shown parameters in the inset. 

by a leveling off of the probability for larger randomness 
(a > 0.75). This is precisely what is observed in Fig. 2 
(due to destruction of band-gap structure) , in Fig. 4 (due 
to creation of semi-localized modes) and Fig. 6 (due to 
both). Finally, we note that an inverse (adjusted) Gaus- 
sian fits well to the obtained data in Fig. 8, as indicated. 
The Gaussian shape of this behavior perhaps underlies 
the observed Gaussian behavior in Fig. 2. 



In this paper we have studied the transmissional prop- 
erties of random two dimensional photonic crystals using 
FDTD method, in two experimentally relevant scenar- 
ios. First we study the effect of randomness in a 90° 
bent wave-guide and secondly we study transmissions 
along the edges of a centrally pulsed photonic crystal. 
By plotting various energy (and field) diagrams as func- 
tions of time and randomness, we have identified var- 
ious interesting effects in the output of such samples. 
We study these properties in a full range of randomness 
and identify two important effects: (i) the creation of 
semi-localized modes and, (ii) the shrinking of photonic 
band-gap. Furthermore, by correlating these two effects 
to overlap probability we have identified the main cause 
of such effects. We suggest that experiments along our 
study be done in order to provide experimental verifica- 
tion of our results. 

Finally, it has been shown that there is a correspon- 
dence in the active and passive modes in random me- 
^jg^27,28^ Therefore, a possible extension of our present 
study is to introduce gain in our PC. Here, one would 
expect that the semi-localized modes studied here would 
also play a crucial role there. This might have important 
implications for the study of random lasers. 
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